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Abstract 

The properties of asymmetric nuclear matter have been investigated in a relativistic Dirac- 
Brueckner-Hartree-Fock framework using the Bonn A potential. The components of the self- 
energies are extracted by projecting on Lorentz invariant amplitudes. Furthermore, the optimal 
representation scheme for the T matrix, the subtracted T matrix representation, is applied and the 
results are compared to those of other representation schemes. Of course, in the limit of symmetric 
nuclear matter our results agree with those found in literature. The binding energy Eb fulfills 
the quadratic dependence on the asymmetry parameter and the symmetry energy is 34 MeV at 
saturation density. Furthermore, a neutron-proton effective mass splitting of m* < m* is found. 
In addition, results are given for the mean-field effective coupling constants. 

PACS numbers: 21.65.+f,21.60.-n,21.30.-x,24.10.Cn 
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I. INTRODUCTION 



Symmetric nuclear matter has been studied extensively. The conventional nonrelativistic 
approach to nuclear matter, the BHF fBrueckner-Hartree-Fock) theory, goes back to earlier 



works by Brueckner and others 



mm 



. A breakthrough was achieved when the 



first relativistic (Dirac-) Brueckner-Hartree-Fock (DBHF) calculations were performed in 

priQ 

the eighties [2|, |8|, |9|] . It could describe remarkably successfully the saturation properties of 
nuclear matter, saturation energy and density of the equation of state (EOS). 

Investigations of asymmetric nuclear matter, however, are rather rare. The breaking 
of isospin-symmetry complicates the problem considerably compared to symmetric or pure 
neutron matter. Some older studies within the nonrelativistic Brueckner scheme can be 
found in Refs. 



IG 



In Ref. PJ 



a calculation in the nonrelativistic Brueckner scheme is 
presented for the Paris potential. Furthermore, the properties of isospin-asymmetric nuclear 
matter have been investigated in the framework of an extended nonrelativistic Brueckner 



approach in Ref. 



Ref. 
model 



14^15^16 



17 



m 



In addition, relativistic Brueckner calculations are performed in 



IS 



2G 



2l|. Finally, we mention two calculations within the a — lj 



The investigation of asymmetric matter is of importance for astrophysical studies such as 



the physics of supernova explosions 24] and of neutron stars j25| , neutron-rich nuclei 



26 



and their collisions 28]. The interest for the structure of neutron- rich nuclei and the colli- 
sions between these nuclei is only of recent date, because data for asymmetric nuclei were 
scarce in the past. However, this situation changes with recent advances in the development 
of high-intensity radioactive beam facilities that will produce nuclei with large neutron ex- 
cess. Hence, systematic theoretical studies of asymmetric nuclear matter are becoming more 
important. 

An important issue is the determination of the precise form of the nucleon self-energy. To 
determine the Lorentz structure and momentum dependence of the self-energy, the positive- 
energy-projected in-medium on-shell T matrix has to be decomposed into Lorentz invariant 
amplitudes. Because of the restriction to positive energy states ambiguities 2^ arise, be- 
cause pseudoscalar (ps) and pseudovector (pv) components can not uniquely be disentangled 
for on-shell scattering. However, with a pseudoscalar vertex the pion couples maximally to 
negative energy states which are not included in the standard Brueckner approach. This 
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coupling to negative energy states is inconsistent with the potentials used and leads to large 
and spurious nuclear self-eriergy contributions from negative energy states ^30] . Hence, it 



was further demonstrated 



3lj that the conventional pv representation used to cure this 



problem fails. The reason is that pseudoscalar admixtures are still present. Finally, new 
and reliable methods, the complete pv representation Q| and the subtracted T matrix rep- 
resentation [s^l, were proposed to remove those spurious contributions from the T matrix 
for symmetric nuclear matter. In contrast, only the conventional pv representation has been 
applied for asymmetric nuclear matter [18|. 

In this work we describe asymmetric nuclear matter at zero temperature in the relativis- 
tic DBHF approach using the Bonn potential and their bare A^A^ matrix elements V jssj ]. 
Furthermore, the optimal representation scheme for the T matrix, the subtracted T matrix 
representation, is applied. A comparison with other representation schemes is made, such 
as the ps representation, the conventional pv representation, and the complete pv represen- 
tation. In addition, the relativistic Pauli operator is used. Compared to symmetric nuclear 
matter the theoretical and numerical effects are larger because protons and neutrons are 
occupying different Fermi spheres. 

The plan of this paper is as follows. The relativistic DBHF is discussed in Sec. [H] Fur- 
thermore, Sec. mil is devoted to the covariant representation of the in-medium T matrix in 
connection with nucleon self-energy components. The results are presented and discussed in 
Sec. IIVI Finally, we end with a summary and a conclusion in Sec. |^ 



II. RELATIVISTIC BRUECKNER APPROACH 



In this section the relativistic Brueckner approach is discussed. First a short description 
of the relativistic Brueckner approach is given for symmetric nuclear matter. Next the 
modifications for the asymmetric case will be treated. 



A. Symmetric Nuclear Matter 

In the relativistic Brueckner approach the nucleon inside nuclear matter may be viewed 
as a dressed particle as a consequence of its two-body interaction with the surrounding 
nucleons. The in-medium interaction of the nucleons is treated in the ladder approximation 
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of the relativistic Bethe-Salpeter (BS) equation 

T = V + i J VQGGT, (1) 

where T denotes the T matrix and V is the bare nucleon-nucleon interaction. The interme- 
diate off-shell nucleons are described by a two-body propagator iGG. The Pauli operator 
Q accounts for the influence of the medium by the Pauli principle and prevents scattering 
to occupied states. The Green's function G which describes the propagation of dressed 
nucleons in the medium fulfills the Dyson equation 

G = Go + Go^G. (2) 

Go denotes the free nucleon propagator while the influence of the surrounding nucleons is 
expressed by the self-energy E. In Brueckner formalism this self-energy is determined by 
summing up the interactions with all the nucleons inside the Fermi sea F in Hartree-Fock 
approximation 

^ = -iJ {Tr[GT] - GT). (3) 

F 

The coupled set of Eqs. dH)-© represents a self-consistency problem and has to be iterated 
until convergence is reached. Below the approximations which are made in the standard 
relativistic Brueckner approach to solve the coupled set of Eqs. dT}-© are discussed. 



1. Self- consistent spinor basis 

Because of the requirement of translational and rotational invariance, hermiticity, parity 
conservation, and time reversal invariance, the most general form of the Lorentz structure 
of the self-energy in the nuclear matter rest frame is 

E{k, kp) = Ss(A;, kp) - 70 So(fc, fcp) + 7 ■ k T.^{k, kp), (4) 

The Ss, So, and Sv components are Lorentz scalar functions which actually depend on the 
Lorentz invariants k'^,k ■ j and j^, where denotes the baryon current. Hence, the Lorentz 
invariants can be expressed in terms of ko, |k| and kp, where kp denotes the Fermi mo- 
mentum. The components of the self-energy are easily determined by taking the respective 
traces 3, 

^s = \trm, So = ^tr[7oS], = [7 ■ k S] . (5) 
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The presence of the medium leads to effective masses and effective momenta of the nu- 
cleons 



m*{k,kF) = M + ^eJ:,{k,kF), k*'' = k^" + ^eW{k, kp). (6) 

The Dirac equation written in terms of these effective masses and momenta has the form 

[^^k*"" -m*{k,kF)]u{k,kF)=0. (7) 

In the following we will work in the quasi-particle approximation, i.e. the imaginary part of 
the self-energy 53mS will be neglected, to simplify the self-consistency scheme. To determine 
the self-energy only positive-energy states are needed in the relativistic Brueckner approach. 
Therefore, the full nucleon propagator can be replaced in Eq. by the Dirac propagator 

Goik, kp) = [-f^k*^' + m*{k, kF)]2m6{k*^ - m*\k, kF))Q{k*^)Q{kF - |k|). (8) 

Here k denotes the momentum of a nucleon inside the Fermi sea in the nuclear matter rest 
frame. This momentum is on mass shell. Due to the B-functions in the propagator only 
positive energy nucleons are allowed in the intermediate scattering states which eliminates 
the divergent contributions from negative energy states. Thus we avoid the delicate prob- 
lem of infinities in the theory which generall y w ill occur if one includes contributions from 



rall y w m oi 

0,yQ. 



negative energy nucleons in the Dirac sea 

By the introduction of the reduced effective mass and kinetic momentum 

K = K/ (1 + ^-(^' ^f)) ' ^*(^' ^f) = ^*(^' ^f)/ (1 + ^v(/^, kp)) , (9) 
the Dirac equation in the nuclear matter rest frame can be rewritten as 

[7^^" - m*{k, kF)]u{k, kp) = 0. (10) 

In general the reduced effective mass is density and momentum dependent. To simplify the 
calculation, however, the "reference spectrum approximation" is applied in the iteration 
procedure, i.e. the effective mass of the nucleon is assumed to be entirely density depen- 
dent (|k| = kp). The method implies that the self-energy itself is only weakly momentum 
dependent. Therefore, at the end of the calculation one has to verify the consistency of the 



assumption ^ S(|k| = kp) with the result of the iteration procedure. 

The solution of the Dirac equation in Eq. (fTUI) provides the in-medium nucleon spinor 



u\{k, kp) 



\ 



E*{k) 



2m*p 



1 

2A|k| 
V B*(k)+n 



XX, 



(11) 







where E*(k) = ^J'k^ + m*p j56|. x\ denotes a two-component Pauli spinor with A = ±|. 
The normalization of the Dirac spinor is thereby chosen as uxik, kp)ux{k, kp) = 1. Since 
the spinor contains the reduced effective mass the matrix elements of the bare nucleon- 
nucleon interaction become density dependent. This density dependence, which is absent in 
nonrelativistic Brueckner calculations, is actually considered as the main reason for the suc- 
cess of the relativistic DBHF approach concerning the description of the nuclear saturation 



2. Covariant T matrix 

We apply the relativistic Thompson equation [37[ to solve the scattering problem of 
two nucleons in the nuclear medium. Therefore the two-particle propagator iGG in the 
BS equation, Eq. ((T)), has to be replaced by the effective Thompson propagator jsTj. The 
Thompson propagator implies that the time-like component of the momentum transfer in 
V and T is set equal to zero which is a natural constraint in the cm. frame, but not a 
covariant one. The Thompson propagator projects the intermediate nucleons onto positive 
energy states and restricts the exchanged energy transfer by 6{k^) to zero. Thus Eq. (0) is 
reduced to a three-dimensional integral equation. In the two-particle center of mass (cm.) 
frame - the natural frame for studying two-particle scattering processes - the Thompson 
equation can be written as 

T(p,q,x)|e.„.. = V^(p,q) (12) 

+ / Vfnk) ""^^ ^^^'"^^ T(kax) 

J (27r)3 ^*^'^E*2(k)2E*(q)-2E*(k) + ze^ 

where q = (q^^ — q2)/2 is the relative three-momentum of the initial state and k and p are 
the relative three-momenta of the intermediate and the final states, respectively. The Pauli 
operator Q depends on the boost three velocity u into the cm. frame. Hence, the T matrix 
depends on the set of parameters x = {kp,m*p, |u|}. 



The Thompson equation ()12|) for the on-shell T matrix (|p| = |q|) can be solved applying 
standard techniques described in detail by Erkelenz js^. Doing so, the positive-energy 
helicity T matrix elements are determined explicitly via the \JMLS >-scheme. In the 
on-shell case only five, for asymmetric nuclear matter six, of the sixteen helicity matrix 



elements are independent which follows from general symmetries j39|. After a partial wave 
projection onto the \JMLS >-states the integral reduces to an one-dimensional integral 
over the relative momentum |k| and Eq. (fT^ decouples into three subsystems of integral 
equations: the uncoupled spin singlet, the uncoupled spin triplet, and the coupled triplet 
states. To achieve the reduction to the one-dimensional integral equations the Pauli operator 
Q has to be replaced by an angle-averaged Pauli operator Q Isl]. Since the Fermi sphere 
is deformed to a Fermi ellipsoid in the two-nucleon cm. frame, Q is evaluated for such a 
Fermi ellipsoid: 

I k I ^ ^mm 

Q ~ ' 'yu\\i\ ~ ^ ^min < |k| < kmax (1^) 

i^max ^ 1^1 



with kmin = ykp — u'^Ep, k^ax = 'j{uEf + kp), and u = |u|. The integral equations are 
solved by the matrix inversions techniques of Haftel and Tabakin (sj. 

The two-nucleon states are two-fermion states. Due to the anti-symmetry of these states 
the total isospin of the two-nucleon system (I = 0, 1) can be restored by the standard 
selection rule 

(-l)L+s+i = _i (14) 

With help of Eqs. (3.28) and (3.32) in IsQ^ the five independent partial wave amplitudes in 
the helicity representation are obtained from the five independent on-shell amplitudes in the 
I J MLS >-representation. The summation over the total angular momentum yields the five 
on-shell plane-wave helicity matrix elements 



< pX'X\T\x)\qX,X2 >= E (^17^) 4aW < A;A;|TJ'^(p,q,a;)|AiA2 > 



J 



Here 9 is the scattering angle between q and p, with |p| = |q|, while A = Ai — A2 and 
a' = X'l — A2. The reduced rotation matrices dj^'^^id) are those defined by Rose 



(15) 
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B. Asymmetric Nuclear Matter 



In symmetric nuclear matter the Fermi momenta of protons and neutrons are equal. 
However, in asymmetric nuclear matter the protons and neutrons are occupying different 
Fermi spheres, leading to different Pauli operators and corresponding neutron and proton 
effective masses and self-energies. This asymmetry is also reflected in the fact that one 
has three different in-medium interactions of the nucleons. They are treated in the ladder 
approximation of the relativistic Bethe-Salpeter equation 

-^np ~ ^p ^ y QnpGnGpT^p ''^ ^ j QpnGpGnT^p , (18) 

and 

-^np ~ ^np ~^ ^ J ^rip QpnGpGnT^p ~^ ^ j ^'^P QripGnGpT^^ , (19) 

where Tij denotes one of the three different T matrices. The three different bare nucleon- 
nucleon interactions are described by one-boson exchange potentials Vij. In the case of 
neutrons and protons having different effective masses, the helicity matrix elements cease to 
be symmetric under the exchange of these particles leading to an additional sixth indepen- 
dent helicity matrix element for the np interaction. These six independent amplitudes can 
be reduced to five, if one assumes an average mass in the np channel for Vnp- In that case, 
V^p*" is related to V^^'^ by the Fierz transformation, which would not be the case otherwise 
due to the unequal effective masses of neutrons and protons. Therefore, Eqs. (fTHj) and (fT^ 
can be reduced to one equation, 

Tjip ^np ~^ ^ j ^npQ npG nG pT^ip. (^^) 



Next, the two-particle propagator iGiGj has to be replaced by the effective propagator. 
In this work it is the Thompson propagator 'sT^. The effective Thompson propagator is 
given by 

g^. = iG^G, = ^ , (21) 

' E* E* ^ -E* -E* + ie ^ ' 
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where is the invariant mass. 

Furthermore, the Pauh operator Q has to be replaced by an angle-averaged Pauli operator 
Q. The angle-averaged relativistic Pauli operators for the nn and pp interactions are identical 
to the one in the symmetric case which is given in Eq. (fT^ . In contrast to the symmetric 
case, the Pauli operator for the np interaction has to be evaluated for Fermi ellipsoids with 
different sizes. The result for the angle-averaged Pauli operator Q^^ for asymmetric matter 
with a neutron excess is 

Q{juEFn - IPFn) |k| < kmin 

l/2[cos(6'p) - cos(6'„)]e(6'„ - Op) for kmin < |k| < kmax (22) 

with kmin = lluEpn " PFn\, kmax = l^uEpn + PFn), 

Q^^}^ V 7|k||u| ; ' 7|k||u| I - ^ (23) 

otherwise 



Q 



np 



and 



arccos ( ^i^^^) far \^^-^^\ < 1 
V 7|k||u| J ■' I 7|k||U| I - 



(24) 



71 otherwise 
However, the central quantities in the model are the neutron self-energy 

T,n = -i J {Tr[GnTnn] - GnTnn) - i j {Tr[GpTnp] - GpTnp), (25) 

Fn Fp 

and the proton self-energy 

T.p = -i J {Tr[GpTpp\ - GpTpp) - i J (TrfG^T^J - G^T^p). (26) 

Fp Fn 

In Brueckner theory the integrations extend over the Fermi sea of the neutron F„ and of 
the proton Fp. Below, the expressions in Eqs. will be specified for the different 

covariant representations. 



III. COVARIANT REPRESENTATION AND THE SELF-ENERGY COMPO- 
NENTS 

In this section we will consider four different representations of the T matrix: pseu- 
doscalar, conventional pseudovector, complete pseudovector, and subtracted T matrix rep- 
resentation. 
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A. Pseudoscalar (ps) Representation 



The nucleon self-energy components are calculated in the nuclear matter rest frame using 
the trace formulas, Eqs. (0). Since we determine the T matrix elements in the two-particle 
cm. frame, a representation with covariant operators and Lorentz invariant amplitudes in 
Dirac space is the most convenient way to Lorentz-transform the T matrix from the two- 
particle cm. frame into the nuclear matter rest frame [8]. On-shell a set of five linearly 
independent covariants is sufficient for such a T matrix representation in symmetric nuclear 
matter. A linearly independent, however, not unique set of five covariants is given by the 
following Fermi covariants 

S = 1 ® 1, V = ® 7^, T = d'^'^ ® A = 757^^ ® 757^, P = 75 ® 75- (27) 

To circumvent the problem of unphysical contributions in the nucleon self-energy one uses 
antisymmetrized amplitudes F/(|p|,6',x) and F/(|p|,7r — 9,x). Thus, the representation of 



the T matrix is given by 

T\\p\,e, x) = T'''''{\p\,e, x) - Ti'-^(|p|, ^, x), (28) 

where the "direct" and "exchange" parts of the T matrix are defined as 

T^'^^^XIPl.^.a;) = ^[Fsi(|p|,e,a;)S + F^(|p|,^,a;)V + F^(|p|,^,a;)T 

+ Fi(|p|,^,a;)A + F^(|p|,^,x)P], (29) 

and 

T''^^%\pie, x) = {-iY^'\ [i^sdpU - 9, x)s + f^(|pU - e, x)v 

+F^(|p|, TT - 9, x)T + FKIpI, TT - 9, x)A + Fl{\p\,7i - 9, x)P] , (30) 



where the interchanged invariants are defined as |41.] S = SS, V = SV, T = ST, A = SA, 
and P = SP with operator S exchanging particles 1 and 2, i.e. Su{1)(ju{2)t- = u{1)t-u{2)„. 
Here p and 9 denote the relative three-momentum and the scattering angle between the 
scattered nucleons in the cm. frame, respectively. In addition, the five Lorentz invariant 
amplitudes F/(|p|,6',x) with i = {S, V, T,A,P} depend also on x = {kp,m*p, |u|}. Taking 
the single nucleon momentum k = (0, 0, k) along the z-axis, we have for the scalar and 
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vector components of the neutron self-energy 



(27r)3 E*^ 



+\ fj' + F^' - 2i^7 - 2^F - ^ri' (32) 



and 



1 f^vr, d^q q . k 



nnl 



^^^^^ = 4 X + ^s'^ ~ ^^v"" - 2i^r - J 

Corresponding expressions as in Eqs. (jHT|) - (jn^ can be obtained for the components of the 
proton self-energy. 

This representation gives a strong momentum dependence in the self-energy in the sym- 
metric case js^. This strong momentum dependence questions, of course, the validity of 
the 'reference spectrum approximation' used in the present self-consistency scheme. Fur- 
thermore, such a strong momentum dependence leads to unphysical results deep inside the 
Fermi sea since the effective mass drops to values which are close to zero. The strong mo- 
mentum dependence of the self-energy was found to originate mainly from the one-pion 
exchange contribution to the self-energy. Therefore, some representations which have a 
better treatment of the one-pion exchange contribution are given below. 



B. Conventional Pseudovector (pv) Representation 

As discussed in Sec. 1111 Al the set of five covariants used to represent the on-shell T matrix 
is not uniquely defined as long as one works exclusively in the subspace of positive energy 
states js^. Obviously, various alternative sets of five linearly independent covariants exist 
such as conventional pv representation. In this representation the pseudoscalar covariant 
P = 75 ® 75 in the T matrix representation in Eq. is replaced by the pseudovector 
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covariant 

m* + m* m* + m* 

with q = pi — P3 and qo = Ei — E3. It leads to identical on-sliell helicity matrix elements as 
the pseudoscalar covariant. 

Using the vv representation of the T matrix as discussed above the nucleon self-energy 
becomes 0] 

-12F^P + 4Ff + " , " (35) 

T A (^*+^*)2 pyJ' ^ ^ 



E,(k) = i£'"^[-4F{?" + F^--2F?"-2Fr 



2m*2 PVJ 4 Jo (27r)3^ v ■ g v 



and 



nn 
A 



2m;2 PvJ^4 7o (27r)3 |k|2E*p ^ ^ s 

The momentum dependence is still pronounced, because the pion contribution is not yet 
treated correctly as pseudovector as discussed in Ref. 



C. Complete Pseudovector {pv) Representation 

Since Fierz transformation mixes the covariants in the conventional pv representation the 
HF nucleon self-energy is not reproduced when the pseudovector pion exchange potential is 
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used for the bare nucleon-nucleon interaction. Hence, to suppress the undesired pseudoscalar 
contribution of the pion to the nucleon self-energy we have to determine a different pv 
representation of the T matrix. To obtain a 'complete' pv representation the identities 



-(T + T) = S + S + P + P 



are needed. Applying the operator identities 
can be rewritten as 



(38) 

V + V = S + S- P- P (39) 
-(jnni) the "symmetrized" ps representation 



T\\p\,e,x) = gl{\p\,e,x)S-gl{\p\,e,x)S + gii\p\,e,x)iA-A) 
+ gU\p\,0,x)P-gl{\p\,e,x)P, 

where the new amplitudes g} are defined as 



(40) 









-2 


-8 





-2\ 




9i 


1 

~ 4 






-6 
-2 


-16 







2 

-2 




9h 







2 


-8 


4 


2 




\9p) 




1° 


6 


-16 





-y 





(41) 



Due to the linear relations between the amplitudes and gj, the two ps representations 
()28|) and ()40|) of the T matrix lead to identical results for the nucleon self-energy. When we 
replace in ()40|) the covariants P, P by the pseudovector covariants PV, PV, respectively, we 
arrive at the complete pv representation j4l| 

T\\p\,9,x) = gU\p\,9,x)S-gl{\p\,9,x)S + g\{\p\,9,x){A-A) 

+ gly{\p\,9,x)PV - g'~{\p\,9,x)PV, (42) 

with gpY{9) and g~^{9) being identical to gp{9) and gl,{9), respectively. As shown in Ref. 



this representation will be able to reproduce the Hartree-Fock results for the nucleon 
self-energy if we use the pseudovector pion exchange potential as bare nucleon-nucleon in- 
teraction. The self-energy components in the complete pv representation of the T matrix 
are given by 



Ss(k) 



mz 



A Jo (277)3 F*' 



m 



9'i 



1 



*2 _ l.*f^n* 

m*2 + m*^ - 2k*f'q* 



(m* + m*) 



PV-' 



(43) 
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and 



^aJo (27r)3^^s ^^?A E*p(m; + m;)2 V^' ^^^^ 



V - ^ "^'"l q-^^ r^"» 9."- fc. < -F^g; 



(27r)3 |k|2E*p^^s 2m, 
1 f^Fp q-k ^ „p „p 



*2 _ ™*2\ 

, — -a^^^- (45) 



However, as already pointed out in j31|, this representation will not reproduce HF nucleon 
self-energy if other meson exchange potentials than pion exchange potentials are used. 



D. Subtracted T Matrix Representation 

The complete pv representation treated in Sect. 1111 HI will fail to reproduce the HF nucleon 
self-energy if other meson exchange potentials are applied than vr- and (77-) mesons as bare 
interaction. Since the influence of the pion is dominantly given by the single-pion exchange, 
it should be reasonable to treat the bare interaction of the t] and especially of the pion 
separately to the rest. After subtraction of the bare interaction of the 77- and vr-meson V^ ,^, 
the remainder is the subtracted T matrix, 

Tsub = T- K,,- (46) 

The ps representation should be more appropriate for the subtracted T matrix because 
then the higher order contributions of the other meson exchange potentials are not treated 
incorrectly as pseudovector. Thus one chooses the complete pv representation for V^^^ and 
the ps representation for the Tsub to get the most favorable representation. 



IV. RESULTS AND DISCUSSION 

In Fig. n] we present the results for the equation of state for various values of the proton 
fraction Yp = Up/ns using the Bonn A potential. The applied representations are the com- 
plete pv and the subtracted T matrix representation. The two extreme cases are symmetric 

14 




FIG. 1: Binding energy as a function of the baryon density for proton fractions Yp ranging from 
to 0.5. The complete pv (top) and subtracted T matrix representation (bottom) are used. 
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nuclear matter {Yp = 0.5) and neutron matter {Yp = 0). The symmetric nuclear matter 
results agree with those of Ref. 3^. The binding energy curves for intermediate values of 
Yp lie between these two extreme curves. In addition to that, the binding energy Ei, shows 
a nearly quadratic dependence on the asymmetry parameter (3 = Yn — Yp in both represen- 
tations. Furthermore, the symmetry energy is defined with help of the binding energy Eh 
as 

\d^Ef,{nB,f3y 



E, 



sym(^B) — 2 



(47) 



/3=0 



Due to the nearly quadratic dependence of the binding energy on the asymmetry parameter, 
the symmetry energy can be equivalently calculated as 

^sym(^B) = Ehins, /5 = 1) - E^iuB, P = 0) . (48) 
In Table H] the calculated symmetry energy is given at saturation density for each repre- 



representation 




E,y^ [MeV] 


conventional pv 


1.41 


30.80 


complete pv 


1.42 


36.74 


subtracted T matrix 


1.39 


34.36 



TABLE I: The symmetry energy at saturation density for various representations. The Bonn A 
potential is used. 



sentation. To compare the results in Table lU some values from literature are given below. 
The symmetry energy is found to be 25 MeV using the Groningen potential 18[, whereas 
using the Bonn C potential j3| it is found to be 28 Mev. In this context it is worth to 
mention that in Ref. j3] the conventional pv representation was used for the decomposition 
of the T-matrix into Lorentz invariants. The result which the authors of j3l find for Bonn 
C is close to the present value obtained for Bonn A adopting the same projection scheme. 
However, the conventional pv representation suffers from on-shell ambiguities which lead to 
large and spurious contributions for the OPE. Thus we will omit this type of representation 
in the following. On the other hand, in Ref. a sixth amplitude was introduced which 
appears only in the np channel for the case of different neutron and proton masses. For the 
two limiting cases, i.e. symmetric matter and pure neutron matter, this amplitude vanishes 
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identically. Hence the results for the symmetry energy are not affected by this additional 
amplitude. As discussed in [l^, the definition of this amplitude is, however, not unique. 
Since we use an averaged neutron-proton mass for the evaluation of the Vnp matrix elements 
as required by the Bonn code, we work, in contrast to Ref. standardly with five ampli- 
tudes. 

A recent liquid drop model calculation gives a result of 32.65 MeV A recent analysis 
of isovector giant dipole resonance (GDR) data within relativistic mean-field (RMF) theory 
set a range of 34 < Egym ^ 36 MeV |43]. Hence, in the complete pv representation the 
result for the symmetry energy in Table |l] is probably too high. The calculation, based 
on the subtracted T matrix projection scheme, yields a symmetry energy consistent with 
the empirical value. The value of 34 MeV is also in remarkable agreement with a recent 
approach to the nuclear many-body problem based on chiral dynamics in combination with 
QCD sum rules [4^. One has, however, to keep in mind that the saturation density ob- 
tained with Bonn A is slightly higher than the empirical value. Bonn A comes nevertheless 
closest towards the empirical values since Bonn B and C v ield smaller densities but at the 
same time significantly too small binding energies 0, I^J. The comparison performed in 



4J| required therefore a rescaling of the DBHF results in order to adjust the densities. The 
symmetry energy does not change dramatically, i.e. at ub = 0.17 fm~^ the values are 33 
MeV (subtracted T matrix) and 34 MeV (complete pv). 

The density dependence of the symmetry energy is shown in Fig. |2l It is similar as, 
e.g. found in the DBHF calculations of [2l|. Although the EOS for symmetric matter is 
relatively soft (K=230 MeV, subtracted T matrix) the isospin dependence, respectively the 
symmetry energy, is stiff, in particular at high densities. As generally found in relativis- 
tic calculations it is significantly stiffer than in non-relativistic BHF approaches [l^. The 
dependence on the representation, i.e. complete pv or subtracted T matrix, is weak. For 
a quantitative comparison with mean-field phenomenology we compare -Egym in Fig- l2l also 
to a phenomenological results obtained recently within the framework of density dependent 
relativistic mean-field theory The asymmetry parameter 04 is thereby varied from 30 
to 38 MeV. At moderate densities the dependence of -Egym is qualitatively similar to the 
relativistic mean-field parameterizations using 04 = 32 — 34 MeV. However, the density 
dependence of -Egym is generally more complex than in RMF theory. In particular at high 
densities -Egym shows a non-linear and more pronounced increase. For this comparison we 
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FIG. 2: Symmetry energy as a function of the baryon density. The complete pv and subtracted T 
matrix representation are used. In addition parameterizations from density dependent relativistic 

n 

mean-field theory ^45i] are shown where the asymmetry parameter 04 is varied from 30 to 38 MeV. 

have chosen the parameter set from j3] which was obtained at equivalent compression mod- 
uh (K=230 MeV) although the authors of j3] favor a slightly higher incompressibility for 
the symmetric case. 

In Fig. 01 the neutron effective mass is plotted as a function of the proton fraction Yp 
for different representations at ub = 0.166 fm~^. As already observed in the symmetric 
case |32i], the magnitude of the self-energy, respectively of the effective mass, depends cru- 
cially on the projection scheme which is used for the T matrix. However, observables like 
the single-particle potential which are determined by the difference of scalar and vector self- 
energy components are much less affected by the different possible projection schemes since 
these effects cancel in leading order. Here we find that the neutron effective mass decreases 
for the ps representation with increasing proton fraction. This behavior is, however, an arte- 
fact of the ps representation which disappears when a proper decomposition of the T matrix 
is used. Mainly due to a relative large scalar amplitude in the direct part of the T matrix 
in the nn channel compared to that in the np channel, the neutron effective mass increases 
for the conventional pv, the complete pv, and the subtracted T matrix representation with 
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FIG. 3: Neutron effective mass as a function of the proton fraction at ub = 0.166 fm ^ for 
different representations of the T matrix. 



increasing proton fraction. In Fig. 0]the neutron effective mass is plotted as a function of the 
baryon density for various values of the proton fraction using the subtracted T matrix 
representation. Of course, a strong density dependence can be observed. In addition, the 
proton fraction influences the neutron effective mass. The upmost curve is the symmetric 
nuclear matter curve , while the lowest curve is the neutron matter curve. The symmetric 
nuclear matter curve is in agreement with the results of Ref . js^ . For intermediate values 
of the proton fraction Yp the curve lies between these two extreme curves. 

A quantity which is sensitive on the momentum dependence of Ss o.t; is the optical po- 
tential which a nucleon feels inside the nuclear medium. The optical potential is given 
by 

[/(|k|, A;°) = S^dkl) - U^^.m + ^-^''"'^^^^^'^'^'^ (49) 

The strength of the isovector nucleon optical potential, i.e., the Lane potential in 
can be extracted from (f/„ — Up)/ [2(3) at saturation density with (3 = Yn — Yp. From a 
large number of nucleon-nucleus scattering experiments at beam energies below about 100 
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FIG. 4: Neutron effective mass as a function of the baryon density for proton fractions Yp ranging 
from to 0.5 using the subtracted T matrix representation. 
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it can be concluded that this potential decreases with energy. The 
isovector nucleon optical potentials calculated for the complete pv and subtracted T matrix 
representation are shown in Fig. |3| In these cases the optical potential in neutron rich matter 
stays roughly constant up to a momentum of 1.5 fm^^, corresponding to a kinetic energy 
of Ej^in ~ 45 MeV, and then decreases strongly with energy. In mean-field approximation, 
i.e. assuming momentum independent self-energy components, the behavior is different 
resulting in a continuously increasing optical isovector potential [57|. In the present case the 
decrease is caused by a pronounced explicit momentum dependence of the scalar and vector 
isovector self-energy components. The optical isovector potential (f/„ — tAil/(2/3) at zero 



momentum is in good agreement with the empirical value of 22 - 34 MeV j^. Hence, the 
empirical isovector potential extracted from proton-nucleus scattering is well reproduced 
by the present calculation. The DBHF model predicts thereby neutron-proton effective 
mass splitting of m* < m*. The same behavior has been found in 0|. These facts stand 
in contradiction with the conclusion drawn in Ref. which are based on nonrelativistic 
considerations. The same holds if one compares to nonrelativistic Brueckner calculations for 
asymmetric nuclear matter which predict an opposite behavior of the neutron-proton mass 
splitting than their relativistic counterparts. In this context one should, however, be aware 
that nonrelativistic approaches determine usually the effective Landau mass. 



1 dU ^ ^ 



M 



+ 



(50) 



which is in general not identical with the Dirac mass Only in the limit of a constant, 
i.e. momentum independent self-energy including a vanishing spatial component Sv, both 
quantities coincide to leading order in the expansion of the relativistic single-particle energy 
if - in addition - the nonrelativistic single-particle potential shows a parabolic momentum 
dependence. However, the relativistic self-energy is in general momentum dependent and 
also the nonrelativistic potential is more complex 0]. 

The mean-field effective coupling constants can be obtained from the scalar and vector 
self-energies. Such coupling functions parameterize the correlations of the T-matrix in a 
handable way. They can be applied to finite nuclei within the framework of relativistic 
density dependent mean-field theory [^. The four channels are: Dirac scalar isoscalar, Dirac 
vector isoscalar, Dirac scalar isovector, and Dirac vector isovector channel. The neutron 
and proton self-energies, respectively, are calculated at their Fermi momentum. Effective 
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FIG. 6: The effective coupling constants in the isoscalar scalar and vector as a function of the 
baryon density for proton fractions Yp ranging from to 0.5 using the complete pv and subtracted 
T matrix representation. 



coupling constants are then given by 

2 n'"^ + n'"p 



9p ^2 



2 K-n; 



(54) 



where n| and are the scalar and vector densities js^. The results for the isoscalar 
and isovector coupling constants are plotted in Figs. IHl and [71 respectively. The representa- 
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FIG. 7: The effective coupling constants in the isovector scalar gg and vector gp as a function of the 
baryon density for proton fractions Yp ranging from to 0.4 using the complete pv and subtracted 
T matrix representation. 



tions used are the complete pv representation and subtracted T matrix representation. As 
already mentioned before, the magnitude of the self-energy depends crucially on the projec- 
tion scheme which is used for the T matrix. As a consequence the coupling constants are 
also depending on the projection scheme. For the complete pv representation the strength of 
the isoscalar coupling constants clearly decreases as the density increases. Applying the sub- 
tracted T matrix representation - which is the more reliable method - the density dependence 
is less pronounced. At low densities, both the scalar and the vector coupling g^^^ show 
a strong decrease with densities but then the vector coupling stays more or less constant. 
Such a behavior is qualitatively similar to relativistic mean-field lagrangians which include 
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non-linear a self-energy terms. When cast into the form of density dependent couplings 
these models have a slightly decreasing scalar and a constant vector coupling. However, 
when applied to finite nuclei the DBHF isoscalar couplings should be renormalized in order 
to shift the too large saturation density (0.18 fm^'^ with Bonn A) to the phenomenological 
value. After such a procedure the isoscalar DBHF self-energy components are in remarkable 
agreement with the results of ^ based on ChPT -|- QCD sum rules. A renormalization of 
the DBHF results can be motivated by higher order corrections of the hole-line expansion, 
such as e.g. contributions from 3-body forces which are known to shift the saturation density 
towards smaller values jl^j ]. 

If DBHF could perfectly be mapped on mean-field phenomenology the isoscalar coupling 
functions should not depend on the isovector parameter. Since the T-matrix is a complex 
object where isovector and isoscalar contributions mix such an idealization can not be ex- 
pected to hold completely. However, the influence of the density on the coupling constants 
is much larger than that of the proton fraction. For the subtracted T matrix representation, 
the isoscalar coupling constants show also some dependence on the proton fraction. How- 
ever, here as well as in the work in Ref. j3| the dependence of the coupling constants on the 
proton fraction is generally weak. For the subtracted T matrix representation the isoscalar 
vector channel shows the strongest variation as a function of the proton fraction Yp which 
lies in this case between 5 and 10 %. 

The strength in the isovector channel is small compared to that in the isoscalar channel, 
especially for the complete pv representation. Furthermore, a strong dependence on the 
proton fraction is observed for the subtracted T matrix representation. The results for the 
isovector channel in case of the subtracted T matrix representation lie closest to the results 
of Ref. ^1 which are based on the Groningen NN interaction. In particular we observe 
the same increase of the scalar isovector coupling at large densities. The variation of the 
isovector couplings as a function of the proton fraction Yp is generally somewhat larger than 
observed in jlj]. 

V. SUMMARY AND CONCLUSION 

We present a calculation of asymmetric nuclear matter in a relativistic Brueckner frame- 
work using the Bonn A potential. The standard treatment, in which the positive-energy- 
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projected on-shell T matrix has to be decomposed into Lorentz invariant amplitudes, is 
applied. Furthermore, the T matrix is represented by a set of five linearly independent 
Lorentz invariants with an average mass for the np interaction. Since the restriction to pos- 
itive energy states causes on-shell ambiguities concerning the pseudoscalar or pseudovector 
nature of the interaction j^^Ql; some freedom in the choice of the representation exists and 
the relativistic nuclear self-energy can not be determined in a unique way. This ambiguity is 
minimized by separating the leading order, i.e. the single-meson exchange, from the full T 
matrix. Therefore, the contributions stemming from the single- vr and-?] exchange are given 
in the complete pv representation. For the other single-meson exchanges and the remaining 
higher order correlations, the ps representation is chosen. 

With this representation scheme, the subtracted T representation scheme js^, we have 
performed calculations for several observables in asymmetric nuclear matter and compared 
them with the results for other representation schemes, e.g. the complete pv representation 
scheme. The binding energy shows the expected quadratic dependence on the asymmetry 
parameter in both representations. The symmetry energy is found to be 34 MeV at sat- 
uration density, whereas in the complete pv representation the value of 37 MeV is found, 
which is too high. With increasing proton fraction the neutron effective mass increases. In 
addition, the density dependence is stronger than the asymmetry dependence on the neutron 
effective mass. Furthermore, a mass splitting of m* < m* is found in both cases. 

In summary, Bonn A predicts a nuclear equation of state which is soft for symmetric 
matter (K=230 MeV) but has a stiff isospin dependence (-Egym = 34 MeV). Both facts are 
in agreement with constraints obtained from heavy ion reactions where kaon production can 

n n 

)e used to constrain the symmetric part [53|| and pion production [5J| and isospin diffusion 

5^ to constrain the symmetry energy . 

The results are also analyzed in terms of mean-field Dirac scalar-vector isoscalar-isovector 
quantities. In general the strength in the isovector channel is small compared to that in the 
isoscalar channel. However, the found strength in the isovector channel is still significant 
for the subtracted T matrix. Furthermore, a strong dependence on the proton fraction is 
observed in that case. 

Therefore, the results above have consequences for nuclei with large neutron excess and 
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for 



neutron stars. 
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